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, • Abstract 

o 

^SJ We analyse in detail the thermodynamics in the canonical and grand 

canonical ensembles of a class of non-asymptotically flat black holes of the 
J> Einstein- (anti) Maxwell-(anti) Dilaton theory in 4D with spherical symme- 

try. We present the first law of thermodynamics, the thermodynamic analy- 
sis of the system through the geometrothermodynamics methods, Weinhold, 
Ruppeiner, Liu-Lu-Luo-Shao and the most common, that made by the spe- 
cific heat. The geometric methods show a curvature scalar identically zero, 
which is incompatible with the results of the analysis made by the non null 
specific heat, which shows that the system is thermodynamically interact- 
ing, does not possess extreme case nor phase transition. We also analyse the 
local and global stability of the thermodynamic system, and obtain a local 
and global stability for the normal case for < 7 < 1 and for other values 
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of 7, an unstable system. The solution where 7 = separates the class of 
locally and globally stable solutions from the unstable ones. 

Pacs numbers: 04.70.-s; 04.20.Jb; 04.70. Dy. 

1 Introduction 

The discovery that a black hole can radiate a temperature characteristic 
of a blackbody spectrum, was made by Hawking. This opened the possibility 
of constructing a semi-classical thermodynamics for systems of solutions of 
black holes arising from general relativity and its recent modifications. 

An interesting class of black holes is that which comes from the minimal 
coupling in 4 dimensions, from the Einstein-Hilbert action with a scalar field 
and Maxwell field, commonly called the Einstein-Maxwell-Dilaton (EMD) 
theory. This theory, imposing spherical symmetry and a static metric, can 
provide non-asymptotically flat (NAF) solutions. These solutions have sev- 
eral characteristics that have not yet been studied with great detail. More- 
over, we also observe that their thermodynamics is not completely estab- 
lished. This class of solutions was first obtained in [1]. Several general- 
izations have been made, as the case with rotation [2], d-dimensional [3], 
multicenter solutions [4], black branes on the linear dilaton background [5] 
and phantom black holes with degenerate horizon [6]. Here, we will widely 
study the thermodynamics of the class of solutions (2.27) and (2.28) of [7]. 

Several other solutions of non asymptotically flat or non asymptotically 
(A)dS black holes have were studied in detail in [8]- [18]. These solutions 
attract much attention due to the possibility of further elucidate the un- 
derstanding of the AdS/CFT correspondence, regarding the linear dilaton 
solutions, arising in a near horizon limit [19], which can lodge holography 
[20] . One of the generalizations of the usual black hole solutions is for a non 
trivial topology, as in [17], this enables us different thermodynamic proper- 
ties of the usual. The Hawking-Page phase transition [ ] is an important 
thermodynamics property for AdS-type black holes, which does not seem to 
occur to these generalizations. Here, we will also do this analysis as a means 
of comparison with other results previously established. 

There are some methods of analysis in black hole thermodynamics theory 
dubbed as geometrical, because they make use of differential geometry to 
determine thermodynamic properties such as: points at which black holes 
become extremal or they pass through a phase transition and thermody- 
namic stability of the system. One of the first methods was proposed by 
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Rao [22], subsequently developed by other authors [ ]. Later, the work of 
Weinhold [2 ] and Ruppeiner [25] were frequently used for the study of the 
black hole thermodynamics. For analysing the thermodynamics of the nei- 
ther asymptotically flat nor (A)dS class of black holes, we use four methods 
called geometries, the Weinhold [24], Ruppeiner [25], Geometrothermody- 
namics (GTD) [26] and that of Liu-Lu-Luo-Shao [27] methods. Regarding 
the latter two methods, they are in a good concordance when different ther- 
modynamic potentials are chosen, like for example the mass and entropy 
representations, this in virtue of the invariance of the formalism by Legendre 
transformations. Finally, we have mention that the results are independent 
of the particular thermodynamic ensemble considered. 

However, we also have to point out that the GTD method can con- 
tain some inconsistencies when compared with the more usual analysis done 
by the specific heat. Recently it has shown that for the cases Reissner- 
Nordstrom-AdS and (phantom case) anti-Reissner- Nordstrom- AdS black holes 
1 28], the GTD method does not reproduce the results obtained by the specific 
heat method. 

For analysing the thermodynamics of the neither asymptotically flat nor 
(A)dS class of black holes, we use four methods called geometries, as well 
as the more usual method that analyses the specific heat of the black hole. 

The paper is organized as follows. In Section 2, we re-obtain a class 
of non asymptotically flat black hole solutions of the EMD theory and its 
thermodynamics variables. We establish the first law of black holes thermo- 
dynamics and calculate its specific heat. In Section 3, we divide the analysis 
of the thermodynamic system through the methods of specific heat, subsec- 
tion 3.1, geometrotermodynamics, subsection 3.2, Weinhold, subsection 3.3, 
Ruppeiner, subsection 3.4, Liu-Lu-Luo-shao, subsection 3.5 and the local 
and global stabilities in the subsection 3.6. The conclusion and perspectives 
are presented in Section 4. 



2 The field equations and the black hole solutions 



For understanding very well the origin of the parameters and structure 
of the solutions discussed here, we will construct the technique of obtaining 
a classes of solutions previously found, as shown in [6]. 

The action of EMD theory is given by: 

S = J dx^^-g [n - 2 rug^V^Vutp + m e 2 ^F^F^ , (2.1) 
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where the first term is the usual Einstein-Hilbert gravitational term, while 
the second and the third are respectively a kinetic term of the scalar field 
(dilaton or phantom) and a coupling term between the scalar and the Max- 
well fields, with a coupling constant A that we assume to be real. The 
coupling constant 771 can take either the value r/i = 1 (dilaton) or 771 = — 1 
(anti- dilaton). The Maxwell-gravity coupling constant 772 can take either 
the value 772 = 1 (Maxwell) or 772 = — 1 (anti-Maxwell) . This action leads to 
the following field equations: 



Dtp = -± mm \e 2X ^F 2 , (2.3) 
iV = 2r h V tl <pV u <p + 2r/ 2 e 2A ^ Q<7^ 2 - F/F v ^j .(2.4) 

From now, we will make use of the same procedure as in [(>]. Let us write 
the static and spherically symmetric line element as 

dS 2 = e 2 ^dt 2 - e 2a ^du 2 - e 2 ^dn 2 . (2.5) 

The metric function a can be changed according to the redefinition of the 
radial coordinate u. Then, we consider the harmonic coordinate condition 

a = 2/3 + 7. (2-6) 

We will also assume that the Maxwell field is purely electric (the purely 
magnetic case may be obtained by electromagnetic duality transformation 
<p — > —ip, F — > e~ 2Xip * F). Integrating (2.2), we obtain 

F 10 ( U ) = ge -2(A V +2/3+ 7 ) (F 2 = _ 2g 2 e -4/3-4A^ ^ (2J) 

with q a real integration constant. Substituting (2.38) into the equations of 
motion, we obtain the equations of second order 

<p" = -mV2\q 2 e 2uJ , (2.8) 

7" = We 2 ", (2.9) 

(3" = - mq 2 e 2 " , (2.10) 

with 

oj = ry-\ip, J = 7 + /3, (2.11) 
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and the constraint equation 



p 2 + 2/?y - 7?i cp 
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2J 2 2ui 

e —mQ e 



(2.12) 



By taking linear combinations of the equations (2.8)-(2.10), this system 
can be integrated and one gets 



00 ' 



.'2 



(p(u) 
Qe 2 " 



-rjiXj{u) + ipiu + ip 



J 



J2J 



where 



A± = (l±??iA 2 ) 



and the integration constants ipo, ip\ G 
The general solution of (2.14) is: 



Q = rj 2 X + q 2 
., a, b G C. 



(2.13) 
(2.14) 
(2.15) 

(2.16) 



oj(u) 



hi 



\Q\a 1 cosh[a(u — no)] 
Q = 0) 



a(u — uq) (a G 



a G 



hi 
hi 
hi 



, Q G R- 

+ , Q G R+) 
Q(«-m )| (a = 0, QGM+) , 
Qa _1 sin[a(u — Uq)] I (a = ia, a, (J G M + ) 



■v/Qa 1 sinh[a(u — uq)] I ( a G 



(2.17) 



(with no a real constant). The general solution of (2.15) reads: 
J(u) = 



ln|6 -1 sinh[6(«-«i)]| (b G R + ) , 
In \u — tii| (6 = 0), 
In lb- 1 sm[b(u - m)) I (6 = i6; 6 G M + ) 



(2.18) 



(ni real constant). In this way we have for A+ ^ 0, the static and spherically 
general solution of the theory given by the action (2.1): 



( dS 2 = e 2 ^dt 2 - e 2a du 2 - e 2 W , 
a(u) = 2J(u) — 7(ii) , 
P(u) = J(u) - 7 («) , 
j(u) = \ + (oj(u) + Ai^in + \(f ) , 
ip(u) = \^ 1 (-r]iXuj(u) + ipiu + ip ) , 
F = —q e 2uj ^du A (it , 



(2.19) 



where, by fixing the spacial infinity at iii = 0, we have six integrations 
constants (q,a,b,uo,<po,(pi), which obey the following constraint equation 
(2.12): 



\+b z 



m<Pi 



(2.20) 
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Now, we will choose one particular class of the solution of that obtained 
in [6]. According to [6], if we fix the solutions of oj(u) and J{u) to be sinh 
functions, with the non-degenerated horizon (m = n = 1) and uq = 0, and 
making the change of radial coordinate e 2bu = /+, where /+ = 1 — (r + /r), 
we obtain the following solution 

dS 2 = r7(r ~ r+ W - , r ° +7 dr 2 - rd+V-^O 2 , (2.21) 
F = TA /lHldr A , e 2 ^ = ( -) ^ , (2.22) 



2r? 2 r V r 

where the parameter 7 = (1 — 7/iA 2 )/(l+ryiA 2 ). This is an exact solution of a 
spherically symmetric, non asymptotically flat (NAF), static and electrically 
charged black hole, with event horizon r + (r+ > 0). The parameters r + and 
r o ( r o ^ 0) are related to the physical mass 3 and the electric charge by: 

(1-7) 

M = y , r+ , (2.23) 



'^•V^' (2 - 24) 

Here, the parameter ro is related to the electric charge. From here, we 
must take the sign (+) in (2.24), which implies in taking the sign (— ) in 
(2.22), for the Maxwell field. This does not mean that we do not treat 
the solution with the sign (— ) from q (q < 0), but rather, we do not write 
explicitly the two possibilities, with q £ M. However, the calculations are 
valid for the two cases, without any loss of generality. We also notice that 
the phantom cases arising from the choice of the function uj(u) in (2.17), 
being sinh, and we impose 772(1 + f?iA 2 ) > 0. Then, the phantom cases 
are reduced to 771 = —772 = — 1 with 1 < 7 < +00 and r/i = 772 = — 1 for 
— 00 < 7 < —1. But, since we need a positive mass for (2.23), we just have 
the unique phantom case 771 = 772 = — 1 for — 00 < 7 < — 1. Considering the 
normal and phantom cases, one gets the following interval 7 G (— 00, +1). 

Now, we are interested in the geometrical analysis representing semi- 
classical gravitational effects of the black hole solutions mentioned before. 
By semi-classical we mean quantize the called matter fields, while leaving 
classical the background gravitational field. Therefore we will work with the 



3 We use the quasi-local mass defined in [29, 30], because it is a non asymptotically flat 
spacetime. 
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semi-classic thermodynamics of black holes, studied first by Hawking [31], 
and further developed by many other authors [32] . 

There are several techniques to derive the Hawking temperature law. For 
example we can mention the Bogoliubov coefficients [33] and the energy- 
momentum tensor methods [34, 32], by the euclidianization of the metric 
[30], the transmission and reflection coefficients [35, 36], the analysis of the 
anomaly term [37], and by the black hole superficial gravity [38]. Since all 
these methods have been proved to be equivalents [ ], then we opt, without 
loss of generality, to calculate the Hawking temperature by the superficial 
gravity method. 

The surface gravity of a black hole is given by [39] : 



5oo 



2 v / -fl'oo5ii 



(2.25) 



where r + is the event horizon radius, and the Hawking temperature is related 
with the surface gravity through the relationship [38, 31] 

T = £ . (2.26) 

Then, for the NAF black hole case (2.21), we get the surface gravity 
(2.25) as 

7 

- r+ (2.27) 



2r^ 



and the Hawking temperature (2.26), in this case is: 



,7 



as 



T=^. (2.28) 
We define the entropy of the black hole as a quarter of the horizon area 

S=-A = - I y/g^Md4 = TrrJ-^rJr 7 . (2.29) 

4 4 J r=r+ 

From (2.22), we can calculate the electric scalar potential on the hori- 
zon. But in this case, we can not directly integrate the component F\q with 
respect to r, because the spacetime is NAF and provides an infinite contri- 
bution when r — > +oo. This problem is known in literature [4], and also 
happens to black holes in (l + 2)-dimensions [40]. We can solve this problem 
by proposing an electric potential able to provide us the Maxwell field in 
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(2.22) and that satisfies the first law of thermodynamics of black holes. We 
then define the following electric potential on the horizon 



A o = -\A^— + «o, (2.30) 
V 2r? 2 r 

where ao is a constant which fixes the potential gauge, from which we derive 
directly the Maxwell field in (2.22). We will check the first law of thermo- 
dynamics for the solution (2.21). Taking the differential of the mass (2.23), 
of the electric charge (2.24) and of the entropy (2.29), we obtain 



dM = ^—^dr + , dg = \f^dr , (2.31) 
4 V 2ry 2 

dS = tt{1 + -i)rlr l ^dr Q + tt(1 - 7)r* +1 'r~~> 'dr + . (2.32) 
In order to to satisfy the first law of black hole thermodynamics, 

dM = TdS + r] 2 A dq , (2.33) 
it is necessary to fix the gauge with 



»» = ^V^' (2 - 34) 



and then, the electric potential (2.30) is written as 



Note that we introduced a compensatory r\i in (2.33), due to the contribution 
of the negative energy density, in the phantom case, of the field of spin 1, F^ u , 
which provides a work with inverted sign in the first law. For studying the 
thermodynamics of a geometric form, it is useful to write the temperature, 
entropy and the electric potential in functions of the mass and the electric 
charge. To do this, from (2.23) and (2.24), we get 



V 1 + 7 1-7 
which, from (2.28), (2.29) and (2.35) provides the temperature, the entropy 



S 



and the electric potential in functions of the mass and the electric charge 
T = W-^.T^ ^ Ml +7)]^ (2>37) 

7T (1 — 7)^ 

5-37 

S = < SiM 1 "V +7 , Si = 712 2 (i- 7 )7-i, (2.38) 

Ml + 7)] 2 

- jVf - /l + 7 \ , 

4, = -A)— , A = 7? 2 z • 2.39 

Q \ 1 -7/ 

Here, we mention that the temperature and the entropy are always positive, 
since we had taken the case q > 0, and the electric potential is always 
negative. 

One can invert (2.38) for writing the mass in terms of the entropy and 
the electric charge 

M(S,q)=q-— (~) . (2.40) 
Then, one can express the specific heat as 

or equivalently by 

Now, we have in hand all the ingredient for studying the thermodynamics 
of this specific system. In the next section we will analyse the specific heat 
and use four geometric methods for studying the thermodynamics of this 
class of NAF black holes. 



3 Thermodynamics of NAF black holes 

3.1 The usual thermodynamics analysis 

Let us start our study of the thermodynamics of the class of NAF black 
hole solutions of EMD theory, using the analysis of the specific heat. The 
specific heat calculated by the use of the mass is directly proportional to 
the entropy, as shown in (2.41). Then, as we have a single horizon r+, there 
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are no extreme case nor phase transition. So, the geometric methods that 
present acceptable results must reproduce the result of the specific heat. In 
general, the thermodynamics system presents thermodynamic interaction, 
since both the specific heat and the temperature of the black hole are non 
null, but do not possess the extreme case, nor phase transition. 

In the next subsections, we will study the system via thermodynamic 
methods called geometries. These methods consist in determining geome- 
try, through a metric for the thermodynamic space of equilibrium states, 
resulting in a curvature scalar of this metric. This scalar can determine 
whether the system possesses phase transition and sometimes if there is an 
extreme case or thermodynamic interaction. We will make the calculations 
by using a mathematical software. 

In subsection 3.6 we will study the local stability by the specific heat 
studied here. In the next subsection we analyse thermodynamic system by 
the Geometrothermodynamics method. 

3.2 The geometrothermodynamics method 

The geometrothermodynamics (GTD) make use of differential geome- 
try as a tool to represent the thermodynamics of physical systems. Let 
us consider the (2n + l)-dimensional space T, which coordinates are repre- 
sented by the thermodynamic potential <E>, the extensive variable E a and 
the intensive variables I a , where a = 1, ...,n. If the space T has a non de- 
generate metric G A b(Z c ), where Z c = {$, E a , I a }, and the so called Gibbs 
1-form G = d& — 5 a bI a dE b , with 5 a b the delta Kronecker; then the structure 
(T, 0, G) is said to be a contact riemannian manifold if A (dO) n ^ is 
satisfied [41]. The space T is known as the thermodynamic phase space. We 
can define a n-dimensional subspace E C T, with extensive coordinates E a , 
by the map ip : E ->■ T, with $ = $(E a ), such that ip* (0) = 0. We call the 
space E the thermodynamic space of the equilibrium states. 

We define the metric for the thermodynamic space of the equilibrium 
states E as [42] 

= { EC S) {^m^) dEadEb ' (3 - 43) 

which is invariant under the Legendre transformations. The main interpre- 
tation of the thermodynamic system, made by this method is realised by 
calculating the curvature scalar of the metric of thermodynamic space of 
the equilibrium states (3.43). With the curvature scalar we can identify 
whether there exists thermodynamic interaction in the system for non null 
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scalar; if there is extreme black hole, when the scale and zero at a specific 
point, and finally if the system undergoes phase transitions. 

For the thermodynamic potential being the entropy (2.38), the metric 
(3.43) has to be given by 

dl as) = ( m om + *aj ) {-8AP dM + W dq ) (3 - 44) 

= 2Sf 7 (l - ^M'^q^+^dM 2 + 2S?7(1 + j)M 2 ^-^q^d^A5) 

The curvature scalar of this metric is identically zero, so there is no ther- 
modynamic interaction according to the interpretation of this method [26]. 
This means that the black hole could not radiate, being compared to an 
extreme black hole. This result is not in accordance with that of the specific 
heat, showing some inconsistency here in this method, for this specific case. 

Making the calculus for the metric of E, using the mass (2.40) as ther- 
modynamic potential, on gets 

2 

x (J-) 11 dq 2 . (3.46) 

The curvature scalar of this metric is also identically zero. Once again the 
incompatibility with the analysis of GTD with that of specific heat has been 
obtained. This was first shown in pathological solutions of phantom black 
holes [42], next in (anti) Reissner- Nordstrom- AdS black holes [28], and also 
shown here. 

3.3 The Weinhold method 

One of the first to formulate a geometric analysis for a thermodynamic 
system was Weinhold. His method, known as the Weinhold method [24], is 
to define a metric of the thermodynamic space of equilibrium states, through 
the mass as thermodynamic potential. This metric is used for the calculation 
of the curvature scalar Rw-, which determine whether the system possesses 
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phase transitions. The Weinhold metric is given by 
,2 d 2 M n , d 2 Af , d 2 M , , 



79 ^ ( S \ (l +7 )g 1-7 (S^H 

2(1 +7) 3^7 / S \ I^r , 2 . _ 

The curvature scalar of this metric is identically null and this shows that 
the system can not be analysed as having or not phase transition. It is 
well known that this method is often inconsistent with the real physical 
thermodynamic system. In the next subsection we study thermodynamics 
through the method of Ruppeiner. 



3.4 The Ruppeiner method 

Ruppeiner [25] also introduced a metric for defining the geometry of the 
thermodynamic space of the equilibrium states, but in this case, with the en- 
tropy as a thermodynamic potential. Again, this metric provides a curvature 
scalar Rr, which shows if the system possesses thermodynamic interaction 
and phase transitions. The Ruppeiner metric is given by 

= S 2 7 (l - 7)M- 1 "V +7 dM 2 - 25i(l +7)(1 - ^)M~^ q 1 dMdq 
-Si 7 (l + 7)M 1 ~ 7 g- 1+7 dg 2 . (3.48) 

The curvature scalar of this metric is identically null. This means that 
there is no thermodynamic interaction, which does not agree with the result 
obtained by the specific heat. It is also known that this method is often 
inconsistent with the thermodynamic description of the physical system. 

In the next subsection we will make the analysis of the thermodynamic 
system by the Liu-Lu-Luo-Shao method. 



3.5 The Liu-Lu-Luo-Shao method 

Recently, Liu-Luo-Lu-Shao [27], in the same order as Weinhold, has 
formulated a metric of the thermodynamic space of the equilibrium states, 
which provides a curvature scalar that determines whether the system pos- 
sesses phase transitions. This metric is a new thermodynamic metric based 
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on the Hessian matrix of several free energy, which in the case of the 
Helmholtz free energy, is given by 



Al* ATAQ^AA A dT Ad^( dA ° 9T \ AQA u_ 9A ° A 2 

LLLS(F) = -dTdS + dA dq = -QgdS + I — - — I dSdq + — dq 



5(1-7)" \SJ 5(1-7) 

x (Jr) i A ° ~ C 1 + f)TiS) dSdq 
2A -!ftql ( S\ T^y , 2 



The curvature scalar of this metric is identically null. As we had seen, 
this result also is not in accordance with that of the specific heat. An 
inconsistency with respect to the specific heat was shown for the topological 
case of black hole in Horava-Lifshitz gravity [43] . 

We then see that for this class of solutions, all the more usual geometric 
methods that use a metric in the thermodynamic space of the equilibrium 
states, were incompatible with the results of the analysis made by the specific 
heat of the black holes. In the next section we will study the thermodynamic 
stability of this class of black hole solutions. 

3.6 The local and global stability 

We now have to study the stability of this thermodynamic system. With 
respect to this, we may split into two parts, the local stability and the global 
one. 

The local stability is easily analysed by the specific heat (2.41) 

C q = (1 ~ 7) 5 . (3.50) 

7 

Here, we see that for the normal case of EMD (r/i =772 = 1), the system 
is always locally stable for < 7 < 1, i.e. C q > 0, and locally unstable 
for — 1 < 7 < 0. But for the phantom case, the system is always locally 
unstable, since 7 G (—00, —1). 

Now we can analyse the global stability of the thermodynamic system. 
First, we define the Gibbs potential in the grand canonical ensemble 

M 

G = M-TS- 7] 2 A q = M(l - TiSi + t^Aq) = . (3.51) 

1-7 
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So, we see that there is no change of the Gibbs potential, unless for the vari- 
ation of the parameter 7. For the normal case, 7 G (—1, 1), the potential is 
always positive and then, the system is globally unstable. But the phantom 
case, 7 G (— 00,— 1), appears as the opposite of the previous one, and the 
system is globally stable. This is not consistent with our previous analysis, 
so that the Gibbs potential, for the grand canonical ensemble, does not seem 
suitable to do this analysis. 

On the other hand, the analysis made by the Helmholtz free energy in 
the canonical ensemble, shows a good agreement with the results already 
obtained in this work. We define the Helmholtz free energy as 

F = M-TS = M(l-TiSx) = -— — . (3.52) 

1-7 

Now the situation seems to have been reversed. In the normal case, we 
have F < for < 7 < 1, which leads to a globally stable system, such as 
the local stability made by the specific heat. However, the phantom case is 
always globally unstable. 

To conclude this section, we discuss the results obtained by specific heat 
and the Helmholtz free energy. 

We see that there is no phase transition when we consider a specific 
solution with a fixed value of 7, for whatever the values of the extensive 
and intensive variables. But there is a transition of the solutions locally 
and globally stable to the unstable ones. This is particularly observed by 
observing the solution where 7 = 0, which separates the class of solutions 
locally and globally stable from that which are unstable. The specific heat 
(3.50) diverges in this case and the Helmholtz free energy (3.52) is identically 
null. This can be seen directly by the metric (2.21) that, for 7 = yields 
the temperature T = l/(4\/2vrg), the entropy S = 4\/2vrMg and the electric 
potential Ao = —M/q. This shows that M = S/(Ay/2irq) = ST, which, 
substituting in (2.41), leads to a divergence and in (3.52) yields zero. Then, 
as C q > and F < for < 7 < 1, C q ->■ 00 and F = for 7 = and 
C q < and F > for 7 G (—00, —1), we get three phases with respect to the 
stability. For < 7 < 1 the thermodynamic of this part of class of locally 
and globally stable solutions, for 7 = the system behaves as a point of 
phase transition, and for 7 G (—00, —1), the system is always locally and 
globally unstable. 

Here we see a great similarity between our analysis and generalized so- 
lutions, for example in [44, 45]. In these papers, as well as in some other, 
we have stable and unstable phases, depending mainly on the parameter of 
the coupling between the dilaton and the matter field, which in our case is 
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that of Maxwell field. There exists a maximum value and one minimum to 
determine the stability of the system. For the parameter A in (2.1) (normal 
case 771 = 1), we have a thermodynamic system globally and locally stable 
when —1 < A < 1 (0 < 7 < 1). For other values of A or for the phantom 
case, we have an unstable system. So here we have no a Hawking-Page phase 
transition, as in [44, 45]. 

4 Conclusion 

We re-obtained a class of non asymptotically flat black hole solutions of 
the EMD theory in 4D with spherical symmetry (2.21). We also established 
a real interval for the parameter 7, such that — 1 < 7 < 1 is for the normal 
case and 7 G (—00, —1) for the phantom case. We calculated the extensive 
and intensives variable for this thermodynamic system. The first law of 
thermodynamics is also established (2.33), by fixing the gauge of the electric 
potential, and calculated the specific heat (2.41), or equivalently (2.42). 

We analysed the thermodynamic system by the geometric methods of 
the geometrothermodynamics, Weinhold, Ruppeiner and that of Liu-Lu- 
Luo-Shao. All analysis by the geometric methods provided an identically 
null curvature scalar for the thermodynamic space of the equilibrium states. 
This is incompatible with the analysis made by the usual specific heat, which 
provides an interacting system without extreme case nor phase transition. 

The local stability analysis done by the specific heat and the global 
stability made by the Helmholtz free energy in the canonical ensemble, shows 
that the normal case is locally and globally stable for — 1 < A < 1 (0<7< 
1). On the other hand, for the other intervals of 7, including all phantom 
cases, the system proves unstable. The solution where 7 = separates the 
class of locally and globally stable solutions from the unstable ones, but 
this can not be a Hawking-Page phase transition. Recently Myung [46] has 
found similar results to ours here, in the case of the Lanczos-Lovelock theory 
in d-dimensional AdS spacetimes. 

We observe that the methods of geometrical analysis of the thermody- 
namic system proved incompatible with the specific heat due to the specific 
form of the potentials and the intensive variables are written in terms of ex- 
tensive variables. With this, our perspective is to change the metrics of the 
thermodynamic space of equilibrium states of these methods, in order to fit 
them with the particularity of this thermodynamic system, as had already 
been suggested for the geometrothermodynamics case in [42, 28]. 
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